Crossover to non-Fermi-liquid spin dynamics in cuprates 
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The antiferromagnetic spin correlation function Sq, the staggered spin susceptibility XQ and the energy 
scale ljfl ~ Sq/xq are studied numerically within the t-J model and the Hubbard model, as relevant to 
cuprates. It is shown that lopl, related to the onset of the non-Fermi-liquid spin response at T > ujfl, is very 
low in the regime below the 'optimum' hole doping c/j < ~ 0.16, while it shows a steep increase in the 
overdoped regime. A quantitative analysis of NMR spin-spin relaxation-rate I/T2G for various cuprates reveals 
a similar behavior, indicating on a sharp, but continuous, crossover between a Fermi-liquid and a non-Fermi- 
liquid behavior as a function of doping. 
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The understanding of the phase diagram of cuprates con- 
tinues to exemplify one of the major theoretical and exper- 
imental challenges |1]. Besides superconductivity (SC) and 
antiferromagnetic (AFM) ordering, several regimes with dis- 
tinct electronic properties have been identified within the nor- 
mal metallic phase. The behavior of spin degrees of free- 
dom, which are the subject of this paper, has been intensively 
studied using the inelastic neutron scattering (INS) jQ|3j] and 
NMR relaxation experiments f?]. They clearly reveal that in 
underdoped cuprates magnetic properties are not following 
the usual Fermi-liquid (PL) scenario within the metallic state 
above the SC transition T > Tc. 

Within a normal FL one expects a dynamical spin suscep- 
tibility Xq(^) to be T- independent at low T,u!. On the 
contrary, INS results show that q-integrated spin suscepti- 
bility exhibits in a broad range of oj and T an anomalous, 
but universal behavior Xl(^) /((jj/r),^rst^ established 
in La2-xSra^Cu04 (LSCO) at low doping |3, 5]. This be- 
havior can be even followed to lowest T in YBaCuaOe+a; 
(YBCO), where Tc has been suppressed by Zn doping 
At the same time, low-energy INS reveals at low T the sat- 
uration of the inverse AFM coiTelation length k = 1/^, at 
■ least in YBCO H and in LSCO at low doping Jill . Anoma- 
' lous T-dependence of ^■^Cu NMR spin-lattice relaxation rate 
1/Ti and of the spin-spin relaxation rate 1 / T2(7 in underdoped 
cuprates is in general compatible with INS |4||, in particular 
l/{TiT) oc x'Li^,T)/uj\^^o oc l/T, in contrast to a T- 
independent value (Korringa law) in a normal FL. 

I On the other hand, cuprates at optimum doping and, more- 
over, in the overdoped regime show a strong reduction of the 
spin response at low energies uj. This is evident from the loss 
of INS intensity in the normal state (as well as in a weak 
resonant peak for T < Tc) and low NMR relaxation rates 
l/Ti, 1/T2G- At the same time, NMR confirms the approach 
to the normal FL behavior, l/(TiT) ^ const, and I/T2G ^ 
const, yi]. There are other indications that the normal FL be- 
havior is approached in the overdoped regime. Recently, the 
angle-resolved photoemission spectroscopy (ARPES) on BiS- 
rCaCuO (BSCCO) system gave evidence for the existence of 
coherent electronic excitations for T < Tx at higher doping 



|7], i.e., the FL-like phase is found in the normal state only in 
the overdoped regime where Tx shows a steep increase with 
hole doping c/j. Intimately related to the onset of the FL-like 
spin response is also the observation that in cuprates doped 
with nonmagnetic Li and Zn the impurity-induced spin sus- 
ceptibility varies as oc 1/(T + Tk), i.e., with a Kondo-like 
behavior with a characteristic temperature Txich) 1 8], where 
Tk in the underdoped regime, whereas it shows a strong 
increase in the overdoped regime. 

From the point of theoretical understanding, an approach 
to a FL behavior in the overdoped regime far from a metal- 
insulator transition seems plausible, nevertheless a solid the- 
oretical evidence is still missing. A crossover from a strange 
metal to a coherent metal phase is, e.g., predicted within the 
slave-boson approach |9]. Frequently invoked interpretation 
is given in terms of the quantum critical point (QCP) at opti- 
mum doping c,* (masked, however, at low T by the SC phase), 
dividing the FL phase at Ch > cj^ and a (singular) non-Fermi- 
liquid (NFL) metal at Ch < ct- While such a concept is well 
established in spin systems |10], its application to metallic 
cuprates is controversial due to the absence of a critical length 
scale ( e.g., ^(T — > 0) — ^ 00). Low energy spin dynamics as 
emerges from INS and in particular from NMR experiments 
has been extensively analysed within the phenomenological 
theory describing a FL close to an AFM instability. In 
the latter approach the spin-fluctuation energy in fact plays the 
characteristic FL scale, as discussed furtheron. 

The present authors recently showed that an anomalous 
Lo/T scaling, as observed at low doping, emerges from a 
gen eral approach to Xq(w) under a few basic requirements 
1I12II : a) collective Q — (tTjTt) AFM mode in the normal 
state is overdamped 7 > o^q, b) equal-time correlations 
Sci ~ (S^qSq) and the corresponding inverse correlation 
length K are finite and saturate at low T. A nontrivial ujq (T) 
dependence then follows from the fluctuation - dissipation re- 
lation, 

(note that we use h — 1 and define Xq(w) in units of g^/i^). 
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leading to a U)/T scaling for T > u)p ^ 76"^^ where 
C oc |fl2|. Within such an approach it is natural that 

cjp > is finite within the whole paramagnetic regime. Nev- 
ertheless, due to strong dependence on parameters, in partic- 
ular on C, ujp{ch) can show quite a sharp crossover from very 
small values in the underdoped regime to a large increase in 
overdoped systems, consistent with experimental indications. 

In this paper we present numerical results for the doping 
dependence of the FL scale cupL within models relevant to 
cuprates, i.e., the planar t-J model and the Hubbard model. 
We furthermore compare these quantities with the ones ex- 
tracted directly from NMR-relaxation and INS experiments 
on cuprates. One possibility is to get ujpL from the full T- 
dependence of various magnetic quantities, in particular from 
static xq{T) and Sq{T). It is evident that in the NFL regime 
T > LOpL a relation follows from Eq. Q, 
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which evolves into the 'classical' relation for A -C S'q. Note 
that A(r) arises from Eq. Q as the integral over the large-w 
tail Xq{^ > T). We are interested in the low-T regime in the 
paramagnetic phase where Sq{T) already saturates. The sat- 
uration is quite evident from the numerical analysis of various 
models Ili llSll . Eq. indicates that even constant Sq can be 
compatible with strongly T-dependent xq{T) which seems 
to be the essence of the NFL regime in cuprates. In contrast, 
one expects a finite xq {T 0) within the FL regime. 

The characteristic energy scale of spin fluctuations is given 
by uipl{T) — Sq/xq{T) with the corresponding T = 
limit uip]:^{0). The latter can be calculated from T ~ nu- 
merical results. Note that ^^^^(O) = (uj) is just the first fre- 
quency moment of the shape function Xq{^,T — 0)/uj for 
u! > 0. On the other hand, we extract ujpL also from experi- 
ments, in particular from NMR 1 /T2G relaxation data, which 
give rather straightforward information on xq(^)- 

Let us first consider the t-J model. 
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with the nearest neighbor hopping on a square lattice, which 
we analyse for J/t ^ 0.3, as relevant for cuprates (for com- 
parison with cuprates we use also t ~ 400 meV). Results for 
Sq (T) and xq (T) are evaluated using the finite-T Lanczos 
method (FTLM) lll4ll . In this way we analyse systems with 

= 18 sites for arbitrary hole doping Ch = N^/N, and with 
Ch < 3/20 for N = 20. It should be also noted that FTLM 
results are rather insensitive to finite-size effects for T > Tfs, 
whereby for systems considered Tfs ^O.lt L14il . 

In Fig. 1 we present results for x — 4Txq as a function of 
Cfi for various T > Tfs- Note that the limiting value within 
the t-J model is x{T 00) = 1 — Ch- Two distinct regimes 
become immediately evident from Fig. I. The crossing of 
curves x{cji) with different T can be used as the definition 
of the 'optimum' doping cjj ~ 0.16, whereby it is indicative 



that the same value is obtained analysing cuprates with high- 
est TejUl. In the underdoped regime x increases by lowering 
T (down to reachable T ^ Tfs) and appears to saturate to the 
NFL behavior, Eq. Q, consistent with the anomalous uj/T 
scaling |12]. On the other hand, at Ch > cj^ the tendency 
of x{T) is opposite. I.e., xq(^) saturates for T < J, indi- 
cating a 'normal' FL behavior. If x{ch) curves would, even 
for lowest T, indeed cross at Ch ~ c^, we would have been 
dealing with a singularity resembling a QCP with diverging 
Xq{T ^ 0) oc 1/T. Moreover XQ(r ^ 0) would be diver- 
gent in the whole regime Ch < c^. Although present results 
cannot exclude this possibility, the deviation visible at lowest 
T = 0.1 < is more in accord with a crossover between FL and 
NFL regimes. 
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Figure 1: (color online) AFM susceptibility 'ITxq vs. doping within 
the t-J model for various temperatures T/t. 

In Fig. 2 we show corresponding FTLM results for 
'-^Fhich) at various T < J. Note that in this regime Sq{T) is 
essentially T-independent, and the values agree very well with 
the T ~ results obtained via the usual Lanczos technique for 
the ground state (g.s.). The latter approach allows the calcula- 
tion of Sq{T = 0) and xq{T — 0) also for somewhat larger 
systems, i.e., for = 20 at all Nh and for iV = 26 at Nh < 2. 
T — results for Sq are shown in the inset of Fig. 2 and over- 
all follow surprisingly well the linear variation 1/Sq = Kch 
with K ^ 15. In contrast to 5q(c;i), the FL scale Wi^^ reveals 
a nonuniform variation with doping. Again, for Ch > c^ojfl 
is already T-independent for T < J, or at least approaching 
finite ll1pl{Q). In the regime Ch < cjj ujpl{T) is decreas- 
ing with T, so that we can establish only an upper bound 
for ojpL- In the same Fig. 2 we plot also results for ujpl{Q), 
evaluated directly via the T — calculation for largest avail- 
able systems. In the overdoped regime the general agreement 
with the FTLM is evident. As expected, in the underdoped 
region obtained ujpl{Q) seem to be consistently smaller that 
ujpl{T > 0) values, whereby a decrease of ujpl{0) with sys- 
tem size is also observed (e.g., values obtained for iV = 26 
systems are smaller than those for N — 18, 20). So we can 



3 



summarize results in Fig. 2 as follows: a) in the overdoped 
regime ujfl ^ Oi{ch — Cho) with Cho ~ 0.12 and a large 
slope a ^ 3.9t ^ 1.6 eV, b) in the underdoped regime our 
results seem to indicate on a smooth crossover to very small 
^FL < J- 
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Figure 2: (color online) FL scale uopL/t vs. c^, obtained for the t-J 
model using the FTLM for T > and the usual Lanczos method for 
T = 0. The inset shows T = results for l/5q vs. Ch- Thin lines 
are guide to the eye only. 



that within the t-J model, in spite of plausibly weaker correla- 
tions within the Hubbard model for L/ = 8 1 1 17]. In Fig. 3 we 
display also the corresponding free fermion result. We notice 
that on approaching the empty band Ce = 1 — ^ both 
curves converge. However, close to half-filling there is a huge 
qualitative difference. 
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Figure 3: (color online) FL scale uopL/t vs. doping c^, as obtained 
via the CPMC method for the Hubbard model with U/t = 8, where 
the dashed line is a guide to the eye. 



The alternative relevant model is the Hubbard model on a 
square lattice, 

^ = E + "■'^■) + C^E^'T^U- (4) 

{ij)s i 

which in the case of strong Coulomb repulsion U ^ t and 
close to half-filling maps onto the t-J model with J = At^ /U . 
We calculate S'q in the g.s. as a function of hole doping Ch 
within the Hubbard model on a square lattice and at [/ = 
8 t using the constrained-path quantum Monte Carlo method 
(CPMC) fl^ . In this method, the ground state wave function 
is projected from a known initial wave function by a branch- 
ing random walk in the overcomplete space of Slater determi- 
nants. Since the method is most efficient in the closed-shell 
cases, we extend our calculations to various tilted square lat- 
tices where the number of sites N is ranging between iV = 34 
to = 164. The susceptibility xq = dniQ/dBQ is calcu- 
lated by computing sublattice magnetization mq induced by 
a small staggered magnetic field Bq . 

Our results for 1 / Sq again reveal a linear variation ^ Kch 
with K ^ 14. Such results are in qualitative agreement with 
previous QMC calculations for U/t = A Uj], where in the 
latter case K ^ 14.3. In Fig. 3 we present corresponding 
'^fl(O). The qualitative behavior of is very similar to the 
result within the t-J model. Fig. 2. In the overdoped regime 
one can again approximate the variation of ujpL as linear, with 
Cho ~ 0.1 and a ~ A.8 t, while for Ch < Cho i^JpL becomes 
very small. Altogether, obtained tJi?/, do not differ much from 



Let us finally estimate xq{T) and consequently lufl di- 
rectly from experiments on cuprates. Within the normal state 
we use the results for the NMR spin-spin relaxation time T2G, 
obtained from the ^^Cu spin-echo decay, related to static Xq 
as 111], 

2G q q 

Assuming that Xq is peaked at commensurate q = Q and 
can be described by a Lorentzian form Xq = XQ'*^/[(q ~ 
Q)^ + K^] (e.g., consistent with INS in YBCO) with k < tt, 
the second term in Eq. Q can be neglected and the form fac- 
tor replaced by -F(Q). This leads to the relation I/T2G ^ 
0.083Ki^(Q)XQ- 1/72G relaxation rates have been mea- 
sured and summarized in Ref. |4!i], i.e., from underdoped to 
optimally doped YBCO with 0.63 < x < 1, underdoped 
YBa2Cu408, nearly optimum doped Tl2Ba2Ca2Cu30io (Tl- 
2223) and the overdoped Tl2Ba2Cu06+5 (Tl-2201), whereby 
the normalization with corresponding F(Q) has been already 
taken into account (see Fig. 8b in Ref. L4J). Note that k rele- 
vant to Xq is the one appropriate for low-w spin dynamics, as 
measured directly by INS (plausibly k < k). For YBCO k{x) 
has been summarized in Ref. |18]. For cuprates considered 
here appropriate hole concentrations Ch have been estimated 
in Ref. 1 15]. Assuming a continuous variation of K{ch) we de- 
termine also K for YBa2Cu408, Tl-2223 and Tl-2201 (for the 
latter we take k = 1.2/ao), not available experimentally. In 
this way, we evaluate xq{T)- Equal-time correlations S'q are 
so far not directly accessible by INS. As shown before they are 
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nearly model independent, so we assume here the t-J model 
results to finally extract corresponding ujfl{T) as presented 
in Fig. 4 for various cuprates. 
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Figure 4: (color online) lopl vs. T, as evaluated from the NMR 
relaxation rate I/T2G |4] and the INS width k for various 
cuprates. The inset shows the extrapolated scales ujpl{Q) and O 
vs. doping Ch- 

For T well above Tc, xq [T] extracted from 1 /T2G follows 
the Curie- Weiss behavior, i.e., xci{T) oc 1/(T + 6). Such 
a behavior emerges also within our analytical approach in the 
scaling regime |12]. Hence, for T > 150K, we can well 
parametrize ujfl{T) = lofl{^){^ + T/Q) and present in the 
inset of Fig. 4 the doping dependence of ll!fl{^) and 8. It 
is evident that both ijJfl{^) and Q reveal a similar behavior, 
which qualitatively and to some extent even quantitatively fol- 
low our result within the t-J and Hubbard models. In partic- 
ular, there is a clear change of scale between the underdoped 
and overdoped cuprates. 

In conclusion we present evidence, based both on numer- 
ical results within t-J and Hubbard models as well as on 
the analyses of NMR and INS experimental data on cuprates, 
that the FL scale lofl exhibits a rather sharp crossover be- 
tween a steep increase in the overdoped regime and very low 
lofl ^ J vci the underdoped regime for Ch < c^. Note that 
in the latter regime within cuprates one can easily reach val- 
ues ojfl (0) smaller than Tc- This can explain why anomalous 
NFL scaling of the spin response as well as of other quanti- 
ties is observed throughout the normal phase at T > Tc- On 
the other hand, the transition to the normal FL is quite abrupt 
in the overdoped regime, at least with respect to the spin re- 
sponse discussed here. 

Our results are well in agreement with other experimental 
evidence for the existence of transition to the FL behavior in 
cuprates. The FL scale Tx, as revealed by recent ARPES 
experiments on BSCCO |7], in particular its doping depen- 
dence Tx{ch) in the overdoped regime, is close to our results 
for ujFL{ch) with an extrapolated Cha ^ 0.1. Similar doping 
dependent scale Tk, analogous to a Kondo scale in metals. 



arises from the analysis of the local-moment susceptibilities 
in YBCO with in-plane nonmagnetic Li and Zn impurities 
1 8]. Experiments show an abrupt and steep increase of Tk 
on approaching the optimum doping. It is plausible that the 
impurity-induced uniform susceptibility is related to local xl 
(and to staggered Xq) in an uniform system, hence Tk seems 
to be related to 0. Needless to say such a relation requires a 
theoretical justification. 

Still, our numerical results cannot exclude the possibility 
of the existence of a QCR From our analysis, the latter can 
be present at the point where ll1fl{^) vanishes on approach- 
ing from the overdoped side, i.e., in our model systems at 
Ch ^ Cho < c,*. An analogous interpretation might follow 
also from experimental values in Fig. 4, as well as from results 
on the Kondo temperature T/f(c/i) |8J. However, the main 
obstacle to such a scenario is that there is no evidence for an 
ordered AFM phase for Ch < Cho, neither from calculated Sq 
within the t-J and Hubbard models nor from experiments. As 
our analysis shows ojf l (0) remains finite throughout the nor- 
mal phase at all dopings Ch down to the onset of the ordered 
AFM phase at c^^*^ < Cho- The experimental distinction be- 
tween the QCP and the present crossover scenario is that in 
principle ujfl{0) > in the normal phase even in the heavily 
underdoped regime, hence one should be able to detect this 
experimentally by suppressing the SC phase, e.g., as investi- 
gated with INS on YBCO system |6]. However, the theory 
II2I1 reveals that ujFLiO) (x uip can be extremely small in the 
underdoped regime. 



[1] M. Imada, A. Fujimori, and Y. Tokura, Rev. Mod. Phys. 70, 
1039 (1998). 

[2] L. R Regnault et at., Physica B 213-4, 48 (1995). 

[3] M. A. Kastner, R. J. Birgeneau, G. Shirane, and Y. Endoh, Rev. 

Mod. Phys. 70, 897 (1998). 
[4] C. Berthier, M. H. Julien, M. Horvatic, and Y. Berthier, J. Phys. 

I France 6, 2205 (1996). 
[5] B. Keimer etai, Phys. Rev Lett. 67, 1930 (1991). 
[6] K. Kakurai et ah, Phys. Rev B 48, 3485 (1993). 
[7] A. Kaminski, etal., Phys. Rev Lett. 90, 207003 (2003). 
[8] J. Bobroff, et al, Phys. Rev Lett. 83, 4381 (1999). 
[9] P A. Lee and N. Nagaosa, Phys. Rev B 46, 5621 (1992). 
[10] S. Sachdev, Quantum phase transitions, (Cambridge University 

Press, Cambridge, 1999). 
[11] A. J. MiUis, H. Monien, and D. Pines, Phys. Rev B 42, 167 

(1990); Y. Zha, V. Barzykin, and D. Pines, Phys. Rev B 54, 

7561 (1996). 

[12] P Prelovsek, I. Sega, and J. Bonca, Phys. Rev Lett. 92, 027002 
(2004). 

[13] R. R. P Singh and R. L. Glenister, Phys. Rev B 46, 11871 
(1992). 

[14] J. Jaklic and P Prelovsek, Adv. Phys. 49, 1 (2000). 

[15] J. R. Cooper and J. W. Loram, J. Phys. I France 6, 2237 (1996). 

[16] S. Zhang, J. Carlson, and J. E. Gubernatis, Phys. Rev Lett. 74, 

3652 (1995); Phys. Rev B 55, 7464 (1997). 
[17] J. Bonca and P Prelovsek, Phys. Rev B 67, 085103 (2003). 
[18] A. V. Balatsky and P Bourges, Phys. Rev Lett. 82, 5337 (1999). 



